SUSSEX preprint SUSSEX-AST 98/3-3, [a.stro-ph/9803244 

Galaxy clusters at 0.3 < z < 0.4 and the value of Qq 



Pedro T. P. Viana 1 ' 2 and Andrew R. Liddle 3 * 

1 Centro de Astrofisica da Universidade do Porto, Rua das Estrelas s/c, J^lOO Porto, Portugal 

2 Instituto Superior da Maia, Castelo da Maia, 4470 Maia, Portugal 

3 Astronomy Centre, University of Sussex, Falmer, Brighton BN1 9QJ, UK 



oo 

OS 

o> 

> 
o 

CO 

> 

CN 
CO 

o 

oo 

0\ 



6 
h 

03 



- ■ 1 February 2008 



ABSTRACT 

The observed evolution of the galaxy cluster X-ray integral temperature distribution 
function between z = 0.05 and z — 0.32 is used in an attempt to constrain the 
value of the density parameter, O , for both open and spatially-flat universes. We 
estimate the overall uncertainty in the determination of both the observed and the 
predicted galaxy cluster X-ray integral temperature distribution functions at z = 0.32 
by carrying out Monte Carlo simulations, where we take into careful consideration all 
the most important sources of possible error. We include the effect of the formation 
epoch on the relation between virial mass and X-ray temperature, improving on the 
assumption that clusters form at the observed redshift which leads to an overestimate 
of Oo. We conclude that at present both the observational data and the theoretical 
modelling carry sufficiently large associated uncertainties to prevent an unambiguous 
determination of Oo. We find that values of Oo around 0.75 are most favoured, with 
O < 0.3 excluded with at least 90 per cent confidence. In particular, the O = 1 
hypothesis is found to be still viable as far as this dataset is concerned. As a by- 
product, we also use the revised data on the abundance o f gala xy clusters at z = 0.05 
to update the constraint on erg given by Viana & Liddle (1996), finding slightly lower 
values than before. 

Key words: galaxies: clusters - cosmology: theory. 
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et al. 1997; Blanchard & Bartlett 1997; Eke et al. 1998 



The nu mber denoity of rieh eluatcra of galaxica at the prcacnt 

epoch has been used to constrain the amplitude of mass 
density fluctuations on a scale of 8/i -1 Mpc (Evrard 1989; 
Henry & Arnaud 1991; White, Efstathiou & Frenk 1993a; 
Viana & Liddle 1996, henceforth VL; Eke, Cole & Frenk 
1996; Kitayama & Suto 1997). This is usually referred to 
as as, where h is the present value of the Hubble parame- 



ichart et al. 1998) 



Me- 



ter, He 



unita of 100 kmr 



Howe 



the derived 



The best method to find clusters of galaxies is through 
their X-ray emission, which is much less prone to projection 
effects than optical identification. Further, the X-ray tem- 
perature of a galaxy cluster is at present the most reliable 
estimator of its virial mass. This can then be used to relate 
the cluster mass function at different redshifts, ca lculated 
for example within the Press-Sche chter framework ( Press & 



value of ctq depends to a great extent on the present mat- 



Schechter 1974; Bond et al. 1991), to the observed cluster 



ter density in the Universe, parameterized by Oo , and more 
weakly on the presence of a cosmological constant, A. The 
cleanest way of breaking this degeneracy is to include infor- 
mation on the change in the number den sity of rich galaxy 
clusters with redshift (Frenk et al. 1990), the use of X-ray 
clusters for th is pu rpose having been proposed by Oukbir 
fc Blanchard ( 1992 ) and sub s equently further investigate d 
( |Hattori fc Matsuzawa 1995| ; |Oukbir fc Blanchard 1997| ). 



Several attem pts have been made recently, with wildly dif- 
fering results (Henry 1997; Fan, Bahcall & Cen 1997; Gross 
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X-ray temperature function. We can therefore compare the 
evolution in the number density of galaxy clusters seen in the 
data with the theoretical expectation for large-scale struc- 
ture formation models, which depends significantly only on 
the assumed values of Oo and Ao = A/3Hq, the latter being 
the contribution of A to the total present energy density in 
the Universe. 

However, the X-ray temperature of a cluster of galaxies 
is not an easily measurable quantity, as compared to its X- 
ray luminosity. A minimum flux is required, so that there is 
a high enough number of photons for the statistical errors in 
the temperature determination to be reasonably small. Be- 
cause of this, although estimates of the present-day cluster 
X-ray temperature function have been available since the 
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early 90's (Edge et al. 1990; Henry fc Arnaud 1991), the 



change in the cluster X-ray temperature function as we look 
further into the past has been much more difficult to de- 



termine . Estimates for the X-ray temperatures of individual 



clusters with redshifts as high as 0.3 have been available for 
some years (e.g. see David et al. 1993), but only with the 
advent of the ASCA satellite has it been possible to measure 
X-ray temperatures for clusters of galaxies around that red- 
shift in a systematic way, and to go to even higher redshifts. 

The evolution of the cluster X-ray luminosity function 
with redshift, though easier to determine, provides much 
weaker constraints on flo and Ao, due to the fact that the 
X-ray luminosity of a galaxy cluster is not expected to be 
a reliable estimator of its virial mass (e.g. Hanami 1993). 
Though it could in principle provide some indication of the 
change of the cluster X-ray temperature function with red- 
shift, t he problem is that not only is there considerable scat- 



ter in t frie present-day cluster X-ray temperature verses lu- 



minosity relation (David et al. 1993; Fabian et al. 1994) 
but it is also not known how the relation may change with 
redshift, though recently it has been argue d that at least 
up to z = 0.4 it does not seem to evolve (|Mushotzky &" 



Scharf 997 



Allen fc Fabian 199 8 
dat, Bl jmchard fc Oukbir 199^ J 



Reichart et al. 1998; Sa- 



The deepest complete X-ray sample of galaxy clus- 
ters pr esently available is the one obtained from the Em 



Henry et al. 1992). This sample is restricted to objects 
with declination larger than —40° and is flux-limited, with 
Fobs > 1-33 x 10 -13 ergcm _2 s _1 , where F ODS is the cluster 
flux in the 0.3 to 3.5 keV band which falls in a 2'. 4 x 2'. 4 
EMSS detect cell. It presently contains 90 clusters of galax- 
ies, after a few misidentifications were recently removed 



(Gioia & Luppino 1994; Nichol et al. 1997). This is the only 



complete galaxy cluster catalogue beyond a redshift of 0.3, 
and as such unique in providing the means to distinguish 
between different possible value s for tip and Ao. However, 
until the recent effort by Henry ( 1997 ), very few X-ray tem- 
peratures were known for those galaxy clusters in the EMSS 
sample with redshifts exceeding 0. 15 (s ee Sadat et al. 1998 
for a recent compilation). Henry (1997) used ASCA to ob- 
serve all galaxy clusters in the EMSS cluster sample with 
0.3 < z < 0.4 and F obs > 2.5 x 10~ 13 erg cm" 2 s" 1 . The 
resulting sub-sample of 10 clusters has a median redshift 
of 0.32, and the data obtained for each cluster, the X-ray 
flux, luminosity and temperature, can be found in Table 1 
of Henry (1997). 

We will use this data together with the present-day (me- 
dian redshift 0.05) cluster X-ray temperature function. We 
work within the extended Press-Schechter formalism pro- 
posed by Lacey & Cole (1993, 1994), which allows an esti- 
mation of the formation times of dark matter halos. We will 
assume the dark matter to be cold, and consider the cases of 
an open universe, where the cosmological constant is zero, 
and a spatially-flat universe, such that Ao = 1 — Oo- 



2 THE THEORETICAL MASS AND 
TEMPERATURE FUNCTIONS 

The usual means by which the mass function of virialized 
structures can be determined analytically is through the 



Press-Schechter approach, which has been found to repro- 
duce well the mass functions recovered from variou s TV-body 



simu lations (Lacey & Cole 1994 



Eke et al. 1996 ; Tormen 



199£). The comoving number density of galaxy clusters in a 
mass interval dM about Al at a redshift z is given by 



n(M,z) dM 



(1) 



da(R,z) 



2 pb ( 

nM~^(R,z) dM 



exp 



2<r 2 (R, 



dM. 



stein Medium. Sensitivity Survey (EMSS) (|Gioia et al. 1990 , & Cole 1994 ; 



where pb is the comoving matter density, cr(R, z) is the dis- 
persion of the linearly-evolved density field smoothed by a 
top-hat window function on the comoving scale R, such that 
R 3 = 3M/47rpb, and S c is the linear density threshold as- 
sociated with the collapse and subsequent virialization of 
the galaxy clusters. This last quantity depends to s ome ex- 
tent on the geometry of the collapsing structures ( Monaco 
1995), but since rich galaxy clusters are relatively rare, it 
would seem to be a fair assumption to ta ke their collapse 
to be close to spherical ( Bernardeau 1994 ). The analytical 
calculation would then yield <5 C = 1.69 in the case of an 
Einstein-de Sitter universe, with a decrease at most of 5 per 
cent when one goes to a universe with sub-critical density, 
as long as Qg > 0.1, wheth er or not a cosmological constant 
is present (Eke et al. 1996). This is supported by the results 
from A^-body simulations, which prefer S c — 1.7±0.2 (Lacey 



Eke et al. 1996 



Tormen 1998] ). As we want to 



be conservative, we will allow for this margin of variation in 
the value of 5 C , assuming it to be equivalent to a 95 per cent 
confidence interval. 

If we have in mind a particular shape for the power 
spectrum of density perturbations, we can further simplify 
equation (|l]) by writing the derivative in terms of a(R,z). 
As we will be assuming all the dark matter to be cold, the 
value of u(R, z) in the vicinity of 8/i _1 Mpc can be obtained 
to a good approximation through, 



a(R, z) = as(z) 



where 



R 



8k- 1 Mpc 



j(R) = (0.3r + 0.2) 



2.92 + lo gl 



R 



Sh- 1 Mpc 



(2) 



(3) 



and r is the usual shape parameter of the cold dark matter 
(CDM) transfer function. 

Note that 7(-R) is independent of z, reflecting the fact 
that the shape of the power spectrum of density pertur- 
bations in the case of CDM models does not change after 
the epoch of matter-radiation equality. We will assume that 
T = 0.230±g; 334 at the 95 per cent confidence level, for which 
a good fit to the observed present shape of the galaxy and 
cluster co rrelation functions in the vicinity of 8/t" 1 Mp c is 



obtained (Peacock & Dodds 1994; Viana & Liddle 1996). It 



should be mentioned however that when the abundance of 
X-ray clusters over a range of temperatures is used, the pre- 
ferred value for F seems to be closer to 0.10 ( Eke et al. 1996 ; 



Eke et al. 1995). Though in this case T is not as well con- 



strained, there is clearly a significant disagreement between 
the values obtained through the two methods. The source of 
this discrepancy needs to be investigated, and may lie ulti- 
mately in a failure of models whose power spectrum can be 
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parameterized through T to reproduce the observed pertur- 
bation spectrum (see e.g. Peacock 1997). 

The quantity eg (2) is related with its present value 
ag(O) via the perturbation growth law 



a 8 (z) = cr 8 (0) 



fl("(*)) 1 

s(n ) 1 + . 



(4) 



where the appropriate formulae for g(Q) and £i(z), depend- 
ing on whether the universe is open or spatially-flat, can 
be found in VL [equations (8) and (10), and (9) and (11), 
respectively] . 

Using expression (E ) we can now calculate the derivative 
appearing in equation (fl), and substituting we end up with 



n(M,z) AM 



2 p b 



2.92 + log, 



7T M 2 3a(R, z 
R 

8k- 1 Mpc 



■(0.3r + 0.2)x 



(5) 



exp 



2a 2 (R,z) 



AM. 



In order to transform this mass function into the cluster 
X-ray temperature function, one needs to relate the virial 
mass of a galaxy cluster to its X-ray temperatu re. Su ch a 
relation has been obtained analytically by Lilje (1992) (see 
also Hanami 1993), a nd confirmed through hydrodynamical 
iV-body simulations ( Navarro et al. 199q ; Bryan & Norman 
199S3 

1/2 fef* < 6 » 



M v oc n 1/2 x 



7 — 



-3/2 



(l + 2m )~ 3/2 {kvTfl'h 



3/2, 



where z m is the redshift of cluster turnaround, and 



X = 



(7) 



(8) 



with £ the ratio between the cluster and background densi- 
ties at turnaround. This quantity was calculated numerically 
in VL, depending only on Q. = fl(z m ) via % ~ £l~ b ^ n \ where 
6(S1) = 0.76 — 0.25 Q, in the case of an open universe, and 
b(Q) = 0.73 — 0.23 Q, for a spatially-flat universe. 

The radii of turnaround and virialization, respectively 
r m and r v , are related through 



ij/2 



2-17/2 



(9) 



in the case when a galaxy cluster is assumed to be an ideal 
virialized system collapsed from a top-hat perturbation. The 
presence of significant substructure during collapse would 
lead to dynamical relaxation thus making the clusters more 
compact, decreasing the ratio r v /r m . 

The proportionality constant in expression Q can be 
obtained either through analytical derivation assuming hy- 
drostatic equilibrium (e.g. see Bryan & Norman 1998), or 
by using results from hydrodynamical TV-body simulations. 
We choose the latter option, as it provides an estimate of 
the un certainties involved, allowing for the possibility of de- 



of White et al. (19931), carried out in an Einstein-de Sit- 



ter universe. They imply that a galaxy cluster with a X- 
ray temperature of 7.5 keV has a virial mass of M v — 



(1.23 ± 0.32) x 10 
tion redshift of z c 



h~ 



M 



0' 



for an estimated virializa- 



0.05 ± 0.05 ( |Metzler fc Evrard 1994| ; 
Navarro et al. 1995| ) . This corresponds to a turnaround red- 
shift of Zm ~ 0.67 ± 0.08, since z c and z m are related by the 
fact that the time of collapse and subsequent virialization, 
t c , is twice the time of maximum expansion, i m , and in an 



Einstein-de Sitter universe 



t=~Ho 1 (l + z) 



-3/2 



(10) 



Putting all this together we are now able to normalize ex- 
pression (^|), 

3/2 



M v 



1 



(1.23 ±0.33) x 10 
31 - 3 / 2 



1/2 -1/2 
X 



<7-J (l 



( 2 ~) 

1.67 \ 3/2 (_k B T 
+ 2 m / V7 



(11) 



.5 keV 



h' 1 Ma 



where the error is 1-sigma. This normalization of the virial 
mass — X-ray temperature relati on ag rees very well with that 
obtained by Bryan & Norman (199£), who used the largest 
set of hydrodynamical iV-body simulations ever assembled. 

We now have the problem that even after the back- 
ground cosmology is chosen, by fixing Qq and Ao, expression 
(til) depends on the redshift of cluster turnaround, z m . As 
this can be determined from the virialization redshift, z c , us- 
ing the fact that t m = 2t c [for expressions of t as a function 
of z in open and spatially-flat universes see VL, equations 
(26) and (27)], we simply need a means to estimate the dis- 
tribution of the redshifts of cluster virialization at each given 
virial mass. The most well justified method which provides 
this was put forward by Lacey & Cole (1993, 19 94), though 
it may slightly underestimate z c (Tormen 1998). 

Lacey & Cole constructed a merging history for dark 
matter halos based on the random walk trajectories tech- 
nique, and derived an analytical expression for the proba- 
bility that a galaxy cluster with present virial mass M would 
have formed at some redshift z, 



p(z) 



p(w{z)) 



dui(z) 
Az 



where 



p(w(z)) = 2w(z) (f 1 — l) erfc 



- 2) exp 



w{z) 
w 2 (z) 



(12) 



(13) 



and 



w(z) = 



5 C (a(M,0)/a(M,z) - 1) 



(14) 



with / the fraction of the cluster mass assembled by red- 
shift z. Independently of background cosmology we will as- 
sume / = 0.75 ± 0.15, as after this mass fraction has been 
assembled it is expected that the X-ray tempe rature of a 
cluster o f galaxies will not change significantly ( Navarro et 



viation; . from hydrostatic equilibrium due for example to al. 1995 ). We will consider the uncertainty in the value of 



bulk gas motions and turbulence (Norman & Bryan 1998) 



The hydrodynamical TV-body simulations we use are those 



/ to roughly correspond to a 95 per cent confidence inter- 
val. Although the expression for the formation probability 
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given above was obtained for a power-law spectrum of mat- 
ter density fluctuations with index n = 0, while at the clus- 
ter scale n is e xpected to be close to —2 ( Henry£taL^992[; 



Hanami 1993; Oukbir, Bartlett & Blanchard 1997; 



Marke- 



vitch l f)9S| ), numerical results show that p(w(z)) depends 
only very weakly on n (Lacey & Cole 1993). 

The present comoving number density of galaxy clusters 
per temperature interval d{k B T) with a mean X-ray tem- 
perature of k B T, that formed at each redshift z, can now be 
calculated using the chain rule and equation (@) with z = 0, 



n(k B T,z)d(k B T) dz 



3 M 



n(M,z = 0)p(z)dMdz. (15) 



The present cluster X-ray temperature function at k B T is 
obtained by integrating this expression from redshift zero 
up to infinity, with the virial mass M obtained through ex- 
pression ( |ri| ) taking into account the assumed k B T and the 
value of the integration variable z — z c . 

The cluster X-ray temperature function at any redshift 
z, for some temperature k B T, can be obtained by taking the 
point of view of someone placed at that redshift, i.e. transfer- 
ring the conditions prevalent at that redshift to the present 
epoch. For example, one takes fio = fi(z), and changes the 
normalization of expression ([ll]) so that for z = one ob- 
tains the virial mass — X-ray temperature relation that ap- 
plies at the redshift of interest. 

In order to compare with the available data, we will 
actually need to calculate the cumulative or integral clus- 
ter X-ray temperature function, N(> k B T,z), i.e. the co- 
moving number density of galaxy clusters with an X-ray 
temperature exceeding k B T at redshift z. This is obtained 
from the differential cluster X-ray temperature function, 
n(k B T, z)d(k B T), by integrating it from the minimum X- 
ray temperature required up to infinity. 



3 COMPARING OBSERVATIONAL DATA 
WITH THEORETICAL EXPECTATIONS 

We use two pieces of data, the integral cluster X-ray tem- 
perature functions at z = 0.05 and z = 0.32, where these 
are the median redshifts of the X-ray cluster samples we are 
considering. 



3.1 The low-redshift data 

The integral cluster X-ray temperature function at z = 0.05 
can be d eterm ined using the dataset presented in Henry & 
Arnaud ( 1991 ). We actually use an updated version of it, 
with more accurate X-ray flux and temperature determina- 
tions, which has been kindly provided to us by Pat Henry. 
Without taking into consideration the temperature measure- 
ment errors, we derived for the number density of galaxy 
clusters at z = 0.05 with X-ray temperature exceeding 6.2 
keV 



N(> 6.2 keV, 0.05) = 1.80 x 10" 7 /i 3 Mpc" 



(16) 



which as expected agrees with the value obtained by Eke et 
al. (1996,1998) and Henry (1997). It is als o in good agree- 
ment with the results of Edge et al. (199C) and Markevitch 
([l998|). 



However, what not many people take into account is 
that the existence of temperature measurement errors makes 
the above value a biased estimator of the real value of 
N(> 6.2 keV, 0.05) in the Universe.^ This can be easily seen 
by constructing a large number of datasets with the same 
clusters in each as in the z = 0.05 dataset, but where the X- 
ray temperature for each cluster is randomly drawn from a 
Gaussian distribution with mean and dispersion as observed 
for that cluster. This procedure simulates the repetition of 
the temperature measurements, assuming the cluster X-ray 
temperatures originally measured are the actual ones. If one 
now determines N(> 6.2 keV, 0.05) for each cluster dataset 
thus obtained, the mean of the distribution turns out to be 
2.12 x 10~ 7 /i 3 Mpc -3 . Hence, the existence of measurement 
errors in the X-ray temperature determinations leads to an 
overestimation of the real value for N(> 6.2 keV, 0.05) if one 
uses the actual measured temperatures. This can be easily 
understood if one remembers that as the X-ray tempera- 
ture goes up, the cluster number density decreases. Given 
that to a first approximation it is as probable for a cluster 
with actual X-ray temperature below 6.2 keV to have a mea- 
sured temperature above that value, as it is for a cluster with 
k B T > 6.2 keV to have a smaller measured temperature, the 
net effect will therefore be an apparent increase in the num- 
ber density of galaxy clusters with X-ray temperature above 
6.2 keV. We expect 1.80 x W~ 7 h 3 Mpc" 3 to be an overesti- 
mation of the real value for iV(> 6.2 keV, 0.05) in the same 
proportion as 2.12 x 10~ 7 ft 3 Mpc -3 exceeds the assumed 
correct value 1.80 x 10 _7 /i 3 Mpc -3 used in constructing the 
artificial datasets above. Therefore the corrected best esti- 
mate for the number density of galaxy clusters at z = 0.05 
with X-ray temperature exceeding 6.2 keV is 

N(> 6.2 keV, 0.05) = 1.53 x 10 ±0 16 x 10~ 7 h 3 Mpc" 3 , (17) 

where the errors represent 1-sigma confidence levels. They 
were obtained through a bootstrap procedure, which allows 
an estimation of the uncertainty associated with the sam- 
pling variance, where we constructed 10 4 cluster datasets by 
randomly selecting, with replacement, fr om th e original list 
of 25 X-ray clusters in Henry & Arnaud ( 1991 ). The number 
of clusters in each sample is drawn from a Poisson distribu- 
tion with mean 25, in order to include the counting error 
in the uncertainty. Each time a cluster is selected its X-ray 
temperature is estimated by randomly drawing from a Gaus- 
sian distribution with the mean and dispersion observed for 
the cluster. In this way the temperature measurement errors 
also contribute to the total uncertainty. 

In Eke et al. (1998), that X-ray temperature measure- 
ment errors lead to an overestimate of the real value for 
N(> k B T, z) is dealt with through a Gaussian smooth- 
ing of the temperature distribution function predicted in 
each (f2o,^o) cosmology, which is then integrated to give 
N(> k B T, z). 

There are several reasons why we chose to concentrate 
on galaxy clusters with X-ray temperature exceeding 6.2 
keV. The first is that N(> 6.2 keV) best represents the mean 
curve going through the observational points for AT(> k B T), 
both at z = 0.05 and z — 0.32, as can be seen in Figure 2 of 



t We are very much indebted to Alain Blanchard for pointing 
this out to us. 
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Henry (1997). Also, the Press-Schechter framework should 
work best on the largest scales, i.e. for the highest masses 
and X-ray temperatures, as in hierarchical cosmologies these 
are the ones for which the density field has evolved least, 
therefore keeping its gaussianity to a greater extent. Re- 
lated to this is the problem of shear, which starts becoming 
an important factor in the collapse of density perturbations 
as the density field develops, leading to deviations from the 
idealized spherical collapse situation. Another reason is that 
in the normalization of the relation between virial mass and 
X-ray temperature for galaxy clusters, we used hydrody- 
namical iV-body simulations which do not take into account 
a possible (pre-)heating of the intracluster medium due to 
starbursts and supernovae in the galaxies. This effect is still 
quite difficult to model realistically, but the few attempts 
that have been made seem to show that it becomes more 
important as the cluster virial mass decreases. For a galaxy 
cluster whose X-ray temperature would otherwise be 5 keV, 
the heating may increas e the cluster X-ray temperature by 



Navarro et 



as much as 15 per cent (Metzler & Evrard 1994 
al. 199$. 



3.2 The high-redshift data 

The comoving number density of galaxy clusters with X-ray 
temperature exceeding 6.2 keV at z = 0.32 can be calculated 
using the EMSS sub-sample of 10 clusters with redshifts be- 
tween 0.3 and 0.4 and fluxes above 2.5 x 10 -13 erg cm -2 s _1 , 



for which Henry ( (L997| ) obtained the mean X-ray tempera- 
tures through ASCA. We used the data in Table 1 of Henry 
( 1997 ), regarding the X-ray flux and temperature for each 
cluster, to estimate the integral cluster X-ray temperature 
function at z = 0.32. We did this both for open and spatially- 
flat universes, using the estimator 



N{> k B T, 



(18) 



where the sum is over all clusters with fcaT) > k B T, and 
Vmax,i is the maximum volume in which cluster i could have 
been detected at the 4<r level in the EMSS within the redshift 
shell under consideration (0.3 to 0.4 in our case). 

The steps which need to be taken in order to ca lculat e 
these volumes are described in detail in Henry et al. (1992). 
They involve the determination of the maximum redshift 
at which each galaxy cluster could have been detected as a 
function of it s obse rved flux, using equations (1) and (2) of 
Henry et al. (1992). In this calculation one has to compen- 



sate for the fact that clusters of galaxies are extended ob- 
jects, and thus some of their flux will be outside the EMSS 
detect cell. One therefore needs to estimate the typical core 
radius, from which most of the flux originates, of th e galaxy 
clusters in the EMSS sub-sample from Henry (1997). In the 
absence of data specific to this sub-sample, we use the ra- 
tio between the extended and detect cell fluxes estimated 



in Henry et al. (1992) for a sample of 4 galaxy clusters ex- 



tended within the EMSS with a mean redshift of 0.29. They 
find this ratio to be equal to 2.10 ± 0.19, where we will as- 
sume the error to be la. This implies a typical core radius 
around 0.15 h~ l Mpc, depending on the chosen values for 
flo and Ao, which we will assume to remain the same for 
0.3 < z < 0.4. 



The detection volume for a given galaxy cluster is 
then obtained by summing, over all limiting fluxes of the 
EMSS (see Table 3 of Henry et al. 1992) starting at 2.5 x 
10~ 13 erg cm -2 s -1 , the volumes lying between a redshift 
of 0.3 and whichever is the lesser of 0.4 and the maximum 
detection redshift for the cluster. The errors affecting the 
calculation of the detection volumes are thus those associ- 
ated with the flux measured for each galaxy cluster and the 
compensation for the extended nature of galaxy clusters. 

As in the lower redshift case, the calculation of N(> 

6.2 keV, 0.32) using the X-ray temperatures measured for 
the galaxy clusters found between z — 0.3 and z = 0.4 
would lead to the overestimation of N(> 6.2 keV, 0.32) 
due to the presence of errors in the X-ray temperature 
determinations. Again, we correct for this by simulating 
the repetition of the temperature measurements. The ra- 
tio between the mean value obtained for N(> 6.2 keV, 0.32) 
from all the simulated datasets, and the value one gets for 
N(> 6.2 keV, 0.32) using the original dataset, provides an 
estimate for the expected ratio between the latter and the 
real value for iV(> 6.2 keV, 0.32) in the Universe. After per- 
forming this correction, the best estimate for the number 
density of galaxy clusters at z = 0.32 with X-ray tempera- 
ture exceeding 6.2 keV becomes 

N(> 6.2 keV, 0.32) = 3.98 f^ (no) x 10~ 8 /i 3 Mpc" 3 , (19) 

where B(Q ) = 0.09 + 0.38 fi - 0.29 fig if the Universe is 
open and B(fi ) = 0.25 + 0.94 fi - 0.78 0§ if the Universe 
is spatially-flat. These fitting functions have an associated 
error of less than one per cent and are valid for fio between 
0.1 and 1. 

The bootstrap method is the simplest way of simulating 
the procedure involved in the extraction of a sample from 
a given distribution. In our case this distribution is simply 
the population of galaxy clusters in the Universe in the red- 
shift bin 0.3 < z < 0.4 with EMSS X-ray fluxes exceeding 
2.5 x 10~ 13 ergcm~ 2 s" 1 . The bootstrap method therefore 
allows the estimation of the dispersion one would expect to 
obtain in the values measured for some quantity related to 
that population, for example N(> 6.2 keV, 0.3 2), if the type 
of sampling that led to the dataset in Henry (1997) was re- 
peated a large number of times across the sky. 

3.3 The method of comparison 

Let us now assume that in our Universe N(> 6.2 keV, 0.32) 
takes some particular overall value. We would then ex- 
pect this value to be the mean of the distribution func- 
tion assembled with the values that would be measured for 
N(> 6.2 keV, 0.32) if th e type of sampling that led to the 
dataset in Henry (1997) was repeated a large number of 
times across the sky. On the other hand, we would expect 
that the shape of this distribution would be that obtained 
through the bootstrap procedure mentioned at the end of 
the previous subsection. We are therefore now in a position 
to ask the following question. If N(> 6.2 keV, 0.32) took 
some overall value in the Universe, how probable would it 
be to measure the value for iV(> 6.2 keV, 0.32) given by the 
dataset in Henry ( 1997 ), after correcting it for the displace- 
ment due to errors in the X-ray temperature measurements? 
We can then attach, for each value of ilo, a probability of 
the value for N(> 6.2 keV, 0.32) given by the dataset in 
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Henry (1997) being actually measured. The exclusion level 



on each value ol fio is obtained simply by subtracting this 
probability from one. 

In summary the following steps were taken, so that an 
exclusion level can be associated with each fio based on the 
X-ray cluster datasets for z = 0.05 and z — 0.32: 

(i) The Universe was assumed to be either open or 
spatially- flat, with fio taking a value between 0.1 and 1. 

(ii) The best estimate for N(> 6. 2 keV 0.32) in the Uni- 
verse, given the dataset in Henry (|1997|) , was calculated. 
This is the result shown in equation (|19[), after correcting 
for the effect of the X-ray temperature measurement errors. 

(iii) A bootstrap procedure analogous to that described 
for the z — 0.05 data was performed in order to deter- 
mine the expected shape for the distribution function of 
JV(> 6.2 keV, 0.32) if th e type of sampling that led to the 
dataset in Henry (1997) was repeated a large number of 



times across the sky. The number of clusters in each sam- 
ple is now drawn from a Poisson distribution with mean 10, 
and the input observational errors (in the typical ratio of 
extended to EMSS detect cell fluxes at z — 0.29, and in 
the ASCA X-ray fluxes and temperatures) are modelled as 
Gaussian distributed. 

(iv) Using the method described in the section 2, the 
theoretically-expected overall value for JV(> 6.2 keV, 0.32) 
given the assumed fio was calculated. The normalization erg 
of the spectrum was fixed by the low-redshift data using 
equation (jl7|). 

(v) The distribution function for JV(> 6.2 keV, 0.32), de- 
termined through the bootstrap procedure, was modified by 
dividing the values obtained for N(> 6.2 keV, 0.32) by their 
mean and multiplying them by the value determined in (iv), 
so that this value becomes the new mean and the relative 
shape of the distribution is maintained. 

(vi) We calculated the probability of obtaining a value as 
high, or as low, as that determined in (ii) , given the distribu- 
tion constructed in (v). The exclusion level on the assumed 
fio equals one minus this probability. 

It should be noted that in step (iv), the calcula- 
tion of the theoretically-expected overall value for N(> 
6.2 keV, 0.32) for each assumed fio has some theoretical un- 
certainty associated with it. This results from the uncertain- 
ties associated with the values one should use for the nor- 
malization of the cluster X-ray temperature to virial mass 
relation, S c , f and T. There is also the uncertainty associ- 
ated with the observed value for N(> 6.2 keV, 0.05). The 
total uncertainty in the theoretically-expected overall value 
for N(> 6.2 keV, 0.32), for each assumed fio, will be calcu- 
lated through Monte Carlo simulations described in section 
4, which will also explain the procedure used to incorporate 
this uncertainty into the calculation of the final exclusion 
level on each assumed fio. 

In general, it is incorrect to calculate the exclusion level 
associated with some theoretically-expected overall value by 
directly using the distribution obtained through a bootstrap 
procedure, and asking how probable it would be to obtain 
such a high (or low) value given that distribution. In prac- 
tice, however, if the bootstrap distribution is symmetric, the 
exclusion level calculated this way, and in the more correct 
way described above in points (v) and (vi), is nearly the 
same. Nevertheless, in the case of the high-redshift clus- 



ter data, the distribution for N(> 6.2 keV, 0.32) obtained 
through the bootstrap procedure is highly asymmetric with 
a long right-sided tail. Given that high fio models tend to 
predict an overall value for N(> 6.2 keV, 0.32) which is 
smaller than the value determined through the dataset in 
Henry (1997), identifying the mean of the bootstrap distri- 



bution with the latter rather than with the former would 
lead to an exaggerated difficulty for higher fio models to 
reproduce the observations, and thus to a higher exclusion 
level. On the contrary, lower fio models, which tend to pre- 
dict the inverse, would benefit from incorrectly calculating 
the exclusion level. In the next section we will estimate this 
effect. 

In the case of the low-redshift data, though the boot- 
strap procedure yields a distribution for N(> 6.2 keV, 0.05) 
which is best characterized as a lognormal, in practice 
the deviation from symmetry is sufficiently small for it 
to be preferable, given the much simpler calculations in- 
volved, to associate the bootstrap uncertainty with the cor- 
rected value for iV(> 6.2 ke V, 0.0 5) determined from the 
dataset in Henry & Arnaud ( 1991 ), rather than with some 
theoretically-expected overall value. 

In Figure 1 we show the binned probability dis- 
tributions obtained through the bootstrap method for 
N(> 6.2 keV, 0.32), in the cases of fio = 1 and 0.3, where 
the Universe is assumed either open or spatially- flat. The 
distributions were altered so that their mean now coincides 
with the theoretically-expected overall value in those cases. 
It also shows the corrected value for N(> 6.2 keV, 0.32) de- 



termined from the dataset in Henry (1993). The peaks in the 
distributions are not an artifact of the chosen binning, but 
correspond to different numbers of galaxy clusters in each 
bootstrap sample having an X-ray temperature in excess of 
our chosen threshold 6.2 keV. 

The calculations we have just described assume that 
there is no scatter in the relation between cluster X-ray tem- 
perature and luminosity. This is not correct, as mentioned 
previously, and can lead to an increase in the estimated 
value for N(> k^,T, z). This incompleteness problem wors- 
ens as the threshold X-ray temperature ksT is lowered, as 
one starts considering clusters with X-ray fluxes dangerously 
close to the flux detection limit. For the same threshold X- 
ray temperature, the problem is also potentially much more 
serious in the case of the z — 0.32 data than for the z = 0.05 
data. The reason is simply that for the same flux detection 
limit, the faintest clusters that can be detected nearby have 
X-ray luminosities (and thus temperatures) which are con- 
siderably smaller than those of the faintest clusters further 
away. 

The effect of the scatter in the X-ray cluster 
temperature-luminosity relation in the calculation of N(> 
6.2 keV, 0.0 5) is negligible, as the z = 0.05 dataset in Henry 
& Arnaud (1991) is claimed to be nearly complete down to 
at least 3 keV. 

In the case of the z = 0.32 data, it is not clear whether 
the presence of scatter in the X-ray cluster temperature- 
luminosity relation may affect the determination of N(> 
6.2 keV, 0.32). Due to the scatter, there is a finite proba- 
bility that some of the 5 EMSS galaxy clusters with X-ray 
flux below 2.5 x 10~ 13 erg cm -2 s _1 , that were found in the 



redshift range from 0.3 to 0.4 jHenry et al. 1992j ), may not 
only have an X-ray temperature in excess of the lowest X-ray 
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Figure 1. Binned probability distributions for the expected value 
of N(> 6.2 keV, 0.32) in the case of Sl = 1 (dotted lines) and 
Co = 0-3 (full lines). They were obtained from 10 4 realizations 
of the bootstrap procedure described in the text. Dotted lines 
correspond to Qo = 1 and full lines to Qo = 0.3. The arrows show 
the corrected valu e for N(> 6.2keV,0.32) calculated using the 
dataset in Henry (1997) for both cases. In the upper panel we 
assume the cosmological constant to be zero, while in the lower 
panel we show the results for a spatially-flat universe. 



temper ature present in the sub-sample of 10 clusters from 
Henry ( |l997| ), 3.8 keV for MS1512.4, but also in excess of 
our chosen threshold temperature 6.2 keV. The minimiza- 
tion of this possibility was in fact another reason for our 
choice of 6.2 keV as the threshold temperature. For exam- 
ple, at 4 keV th e effect is already expected to be significant 
(|Ekc et al. 199c|). 



We calculated the expected increase in the corrected 
value of iV(> 6.2 keV, 0.32) given the dataset in Henry 
(1997), as a result of the existence of the 5 EMSS galaxy 
clusters mentioned above, by doing 1000 Monte Carlo simu- 
lations where the X-ray temperatures for those clusters were 
estimated via the X-ray cluster t emper ature-luminosity re- 
lation determined in Eke et al. (1998) using the m ore re - 
cent data for the galaxy clusters in Henry & Arnaud (1991). 



The X-ray temperatures for those 5 clusters were obtained 
through the power-law relation 



k B T - 



1/6 



(20) 



where in each Monte Carlo simulation the values of the 
parameters a and b were drawn from Gaussian distribu- 
tions with respectively mean 2.53 and dispersion 0.69, and 
mean 3.54 and dispersion 0.47. The X-ray luminosities in 



the 0.3 to 3.5 keV band for the 5 clusters in question in 
units of 10 44 h~ 2 erg s 1 , L44 , we re de termined using the X- 
ray fluxes given in Henry et al. (1992), with a K-correction 
of 15 pe r cent included independently of the assumed cos- 
mology ( Henry et al. 1992j ; Henry 1997). The maximum 
volumes within which the 5 clusters could have been ob- 
served were calculated in the same way as those for the 
10 higher flu x clu sters, using the EMSS data provided in 
Henry et al. (1992). The final exclusion level for each value 



of fio was calculated simply by taking the mean of the exclu- 
sion levels associated with the new higher estimated values 
for N(> 6.2 keV, 0.32) in the Universe, resulting from each 
Monte Carlo simulation. 

By using the X-ray cluster temperature-luminosity re- 
lation observed for z — 0.05 to estimate X-ray temperatures 
for galaxy clusters where z is within 0.3 to 0.4, we are im- 
plicitly assuming that this relation does not evolve much 
between a redshift of 0.4 and the present. This is supp orted 
by the recent a nalyse s of Mushotzky & Scharf (L997) and 
Allen & Fabian (1998) (see also Sadat, Blanchar d and Ouk- 
bir 1998). However, the actual dataset in Henry (1997) does 
not support this assumption, as a chi-square fit of a power- 
law to the data prefers the situation where the X-ray lumi- 
nosity is practically independent of the X-ray temperature. 
This is driven by the relatively low X-ray temperature mea- 
sured for the galaxy cluster MS2137.3, which, although it is 
by far the brightest cluster in the dataset, has only the 7th 
highest X-ray temperature. Removing this cluster from the 
dataset, the best-fit power-law for the X-ray temperature- 
luminosity relation becomes compatible with the z = 0.05 
one. It was this very strong dependence of the best-fit X-ray 
cluster temperature-luminosity relation on the inclusion or 
not of a single galaxy cluster which led us to choose not to 
estimate the X-ray temperatures for the 5 clusters with the 
lowest fluxes using the data for the 10 clusters that make up 
the dataset in Henry (1997). 

In the end, we found that allowing for the presence of 
scatter in the X-ray cluster temperature-luminosity relation 
when calculating iV(> 6.2 keV, 0.32) has only a small effect, 
at the few per cent level, on the exclusion levels obtained 
for different fio, and does not alter our conclusions. 



4 RESULTS 

The calculation of the theoretically-predicted overall value 
for iV(> 6.2 keV, 0.32) was performed as described in Sec- 
tion 2, with ag obtained from the observed value for 7V(> 
6.2 keV, 0.05), and the uncertainty calculated via a Monte 
Carlo procedure. 

Our first result is the value required for as, as a func- 
tion of both Oo and Ao, so that the observed value for 
N(> 6.2 keV, 0.05) is reproduced. This supersedes the re- 
sult obtained in VL. We find that the best-fitting value is 
given by 



0.56 fig 
0.56 fin 



Open 
Flat. 



This is accurate within 3 per cent for Qo between 0.1 and 1. 

The most important reason why this value is smaller 
than that quoted in VL is the decrease in the assumed num- 
ber density of galaxy clusters at z = 0.05. This results from 
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the revision of the Henry & Arnaud dataset and from the 
correction due to the existence of X-ray temperature mea- 
surement errors, which had not been taken into consider- 
ation in VL. Also, the cluster X-ra y temperature function 
obtained in Henry & Arnaud (1991) h ad been slightly over- 
estimated due to a calculational error (Eke et al. 1996). 



The overall uncertainty in the value of as was calcu- 
lated in the same way as in VL, through a Monte Carlo 
procedure where the sources of error, namely the normal- 
ization of the cluster X-ray temperature to virial mass re- 
lation, the value of 5 C and the value of /, are modelled as 
being Gaussian distributed, and T and N (> 6.2 keV, 0.05) 
as having a lognormal distribution. As in VL, we find that 
for each S7o between 0.1 and 1 the distribution of as can be 
approximated by a lognormal. For open models, the 95 per 
cent confidence limits are roughly given by +20fi Sl ° 



per cent and — 18fi ' 
we have -^Ofi^ 2 ' 0610 " 



.10 p er cem - 
per cent and 



while for flat models 



-180 



o.2io glo n 





The calculation of the uncertainty in the theoretically- 
predicted overall value for N(> 6.2 keV, 0.32), for each as- 
sumed fio, was made using the Monte Carlo simulations per- 
formed in order to calculate the uncertainty in the value of 
as, described above. 

We find that for fio between 0.1 and 1, the distribution 
of N(> 6.2 keV, 0.32) is close to lognormal and, with an 
associated error of less than 4 per cent, its mean is fitted by 



N(> 6.2 keV, 0.32) 



2.67fi ( 7 G(f2o) x 10~ 8 h 3 Mpc" 3 . 



(21) 



where G(fi ) = 0.57 + 2.69 fio - 1.87 fig if th e Universe is 
assumed open and G(fio) = 0.45 + 2.56 fio — 2.12 fig if it is 
assumed spatially-flat. The 95 per cent confidence intervals 
are to a fair approximation given by +170 0° 



0.17+0.31f2 



per 



cent and —62 fi ' ' per cent in the open case, and 
+ 170fi0 ° 6+0 - 42n ° Per cent and -64 fi« 03+016n ° 



per cent m 
the flat case. 

As described in subsection 3.3, we can now build the 
distribution function one would expect to recover if fio took 
a certain value in the Universe and N(> 6.2 keV, 0.32) was 
measured a large number of times across the sky under the 
s ame type of sampling that led to the dataset in Henry 
(1997). Its mean is the theoretically-expected overall value 
for N(> 6.2 keV, 0.32) given in expression ( pl|) for the fio 
under consideration, and the shape of the distribution is 
that obtained through the bootstrap procedure described in 
subsections 3.2 and 3.3. The exclusion level on the assumed 
fio is then given by one minus the probability of measur- 
ing a value for iV(> 6.2 keV, 0.32) as high (or as low) as 
that implied by the dataset in Henry (1997), calculated in 
expression (Jigj) , given the expected distribution. 

However, due to the uncertainties in the estima- 
tion of the theoretically-expected overall value for N(> 
6.2 keV, 0.32), the actual calculation of the exclusion level 
for each fio is not as simple. So that we can obtain it, we need 
to integrate over all possible values for the theoretically- 
expected overall N(> 6.2 keV, 0.32), which we will denote 
u. The overall exclusion is the product of the probability, 
P(u, fio), of each u being the correct overall value one would 
expect for N(> 6.2 keV, 0.32) in the Universe (given the 
assumed fio), and the exclusion level Ex(u) calculated as 
described above for each assumed u, i.e. 
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Figure 2. The absolute exclusion levels for different values of 
Qq in both the open and spatially-flat cases, when the threshold 
X-ray temperature of 6.2 keV is used. 



Exclusion probability of fio 



P(u,n )Ex(u)du (22) 



As mentioned above, the P(u, fio) are lognormal distribu- 
tions with mean given by expression (|2~t]). The dispersion 
can be calculated from the 95 per cent confidence limits. 

In Figure 2 we show the exclusion levels for fio obtained 
in this way. Even for the values of fio for which it is easi- 
est to reproduce the observations, from 0.7 to about 0.8, 
the exclusion level is quite high, around 70 per cent. The 
reason lies with the large uncertainty in the theoretically- 
expected overall value for N(> 6.2 keV, 0.32). Because of it, 
most theoretically-expected overall values end up far away 
from the value for N(> 6.2 keV, 0.32) w hich i s expected in 
the Universe given the dataset in Henry ( 1997 ). A large un- 
certainty in the theoretical prediction is clearly no basis to 
discard models. However, for the high and low fio we are 
aiming to constrain, this effect becomes much less impor- 
tant; the high exclusion levels are caused by most of the 
distribution for the theoretically-expected overall values for 
N(> 6.2 keV, 0.32) being higher (for low fio), or lower (for 
high fio), than the observations. Note that the exclusion lev- 
els are absolute, and not relative as one would obtain from 
the calculation of a likelihood function. 

In the calculation of the theoretically-expected overall 
value for N(> 6.2 keV, z — 0.32), the parameter / repre- 
sents the fraction of the cluster mass at the redshift of clus- 
ter observation, z bs, that had been assembled by the time 
the cluster virialized, z c . We considered this parameter to 
be equal to 0.75, though we allowed for the possibility that 
it could be as low as 0.60 or as high as 0.90, corresponding 
this to something like a 95 per cent confidence interval. We 
modelled this uncertainty by assuming in the Monte Carlo 
simulations that / was Gaussian distributed with mean 0.75 
and dispersion 0.075. 

In order to estimate the effect of changing the assumed 
value for / in the determination of N(> 6.2 keV, z — 0.32), 
we also performed calculations where we treated / as having 
no associated uncertainty. We considered the cases where / 
was equal either to 0.60, 0.75 or 0.90. We found that chang- 
ing / from 0.75 to 0.60 decreased the estimated value for 
N(> 6.2 keV, z — 0.32) by about 10 per cent, while chang- 
ing / from 0.75 to 0.90 increased N(> 6.2 keV, z = 0.32) by 



© 0000 RAS, MNRAS 000, 000-000 



Galaxy clusters at 0.3 < z < 0.4 and the value of f2 9 



from introducing the possibility that z c > z b s ■ As the num- 
ber density of virialized objects evolves faster for the same 
relative change in the value of the dispersion of the density 
field the smaller this value is, the decrease in the required 
value for og, has the effect of enhancing the decrease in the 
value of iV(> ksT, z) as z gets larger. This effect turns out 
to be more important than the expected increase in the value 
of N(> fceT, z) due to higher cluster X-ray temperatures at 
fixed cluster mass resulting from the possibility of z c > z obs . 



5 DISCUSSION 



Figure 3. The expected redshift evolution of 2V(> 6.2 keV, z) for 
Co = 1 ana 0.3. The solid lines show the result obtained using the 
Lacey & Cole method for estimating z c , and the dashed ones the 
result obtained assuming that z c = z bs- Each curve is normalized 
to reproduce the observed value for N(> 6.2 keV, 0.05). Note that 
the divergence at high z is caused by this renormalization; the 
absolute correction is largest at the lowest redshift, where Q(z) is 
smallest. 



around 11 per cent. The impact of removing the uncertainty 
in the value of / from the estimation of the overall uncer- 
tainty in N(> 6.2 keV, z = 0.32) is therefore negligible. 

In all previous uses of the Press-Schechter framework to 
calculate the evolution of the number density of rich galaxy 
clus ters with redshift ( Oukbir fc Blanchard 1992|; |Ekc ct al 



1996j; lOlukbir fc Blanchard 1997t|Henrv 19971 ; |Eke et al. 19981 

Markevitch 1995; Reichart et al. 1998), it has been assumed 
that the redshift of cluster virialization, z c , coincides with 
that at which the galaxy cluster is observed, z b s . In Figure 3 
we compare the value of N(> 6.2 keV, z) obtained using the 
Lacey & Cole (1993, 1994) prescription for the estimation of 
z c with the result of the assumption that z c = z b s - We al- 
ways require that the observed value for N(> 6.2 keV, 0.05) 
is recovered. 

As expected, the difference in the theoretically- 
predicted overall value of JV(> 6.2 keV, z) resulting from 
the two distinct assumptions regarding z c becomes larger for 
Oo = 0.3, reflecting the fact that as Hq goes down galaxy 
clusters tend to form increasingly at an earlier epoch than 
that at which they are observed. We find that neglecting 
the fact that some clusters of galaxies virialize prior to the 
epoch at which they are observed leads to an underesti- 
mation of the predicted degree of evolution in the value of 
N(> ksT, z) for z > z noim , where z nolm is the redshift at 
which N(> ksT,z) is normalized through observations, e.g 
in our case z norm = 0.05. Taking into account the possibility 
that z c may be larger than z a hs therefore requires lower val- 
ues for Qq in order for the high-redshift data on JV(> ksT, z) 
to be reproduced. 

Allowing for z c > z Q b s means that some galaxy clusters 
that otherwise would not be massive enough to reach a given 
threshold temperature k^T can now be counted when calcu- 
lating iV(> ksT, z). In principle this would have the effect 
of increasing the expected value of N(> ksT, z) for any z. 
However, at the normalization redshift 0.05 the higher value 
for iV(> fceT, 0.05) means that a less well developed density 
field at z = 0.05 is required, i.e. a lower value of ag results 



From the above analysis, we conclude that at present it is not 
possible to reliably exclude any interesting value for Slo on 
the basis of X-ray cluster number density evolution alone, 
due to the limited statistical significance of the available 
observational data and to uncertainties in the theoretical 
modelling of cluster formation and evolution. However, we 
do find that values of £Iq below 0.3 are excluded at least at 
the 90 per cent confidence level. Values of fio between 0.7 
to 0.8 are those most favoured, though not strongly. These 
results are basically independent of the presence or not of a 
cosmological constant. 

Our conclusio ns su pport those of Colafrancesco, Maz- 
zotta & Vittorio (1997), who tried to estimate the uncer- 



tainty involved in the estimation of the cluster X-ray tem- 
perature distribution function at different redshifts based 
on its present-day value. They found this uncertainty, given 
the still relatively poor quality of the data, to be sufficiently 
large to preclude the imposition of reliable limits on the 
value of fio- 



O ur res ults disagree with those of Henry ( 1997 ) and Eke 
et al. ( 1998 ) , as they found the preferred Qo to lie between 
0.4 to 0.5, with the Qo = 1 hypothesis strongly excluded. 
This disagreement is mainly the consequence of our focus 
on the threshold X-ray temperature of 6.2 keV, while they 
draw their conclusions based on the analysis of the results 
obtained for several threshold X-ray temperatures. Further 
below we will repeat our calculations assuming a threshold 
X-ray temperature of 4.8 keV, and we will find that when 
we calculate the joint probability of some value for fio being 
excluded on the basis of the results concerning either one or 
both threshold X-ray temperatures of 6.2 keV and 4.8 keV, 
the favoured value for Oo decreases to around 0.55. Some of 
the reasons for our choice of deriving the conclusions solely 
based on the results obtained for the 6.2 keV threshold were 
mentioned at the end of subsection 3.1 and others will be 
detailed below. 

Other less important contributions to the dif ferenc e be- 
tween our r esults and those presented by Henry (1997) and 
Eke et al. (1998) are the different assumed normalization 
for the virial mass to X-ray temperature relation, and the 
corrections in the expected values in the Universe for both 
N(> 6.2 keV, 0.05) and N{> 6.2 keV, 0.32) due to the un- 
certainties in the X-ray cluster temperature measurements. 
Note that changing the mean of the bootstrap distribution 
obtained for N(> 6.2 keV, 0.32) to its theoretically-expected 
overall value in some Qq universe and then calculating the 
exclusion level on the estimated value for N(> 6.2 ke V, 0.32) 
in the Universe given the dataset in Henry (1997), rather 
than just using the original bootstrap distribution to im- 
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Figure 4. The absolute exclusion levels for different values of 
Qq in both the open and spatially-flat cases, when the threshold 
X-ray temperature of 4.8 keV is used. 



pose an exclusion level on the theoretically-expected overall 
value for N(> 6.2 keV, 0.32) in that fio universe, does not 
seem to make much difference. This is a reflection of the 
fact that the bootstrap distributions recovered do not have 
a strongly asymmetric shape. 



Our disagreement with Eke et al. (1998) on the level of 
exclusion of the fio = 1 hypothesis is also due to our much 
larger assumed uncertainty in the theoretically-expected 
overall value for iV(> 6.2 keV, 0.32). 

For the fio = 1 hypothesis to be favoured, one requires 
the lowest possible observed value for N(> 6.2 keV, 0.32). 
This is best achieved if, for the sample of 10 galaxy clusters 
used in its calculation, the X-ray temperatures turn out to 
be on average lower than the assumed mean, and the X-ray 
fluxes higher. A higher ratio between the extended and de- 
tect cell fluxes for the EMSS at z — 0.32 would also help. 
On the theoretical side, the higher one decides the expected 
value for N(> 6.2 keV, 0.32) is, the more compatible with 
the data the fio = 1 hypothesis becomes. This can be best 
achieved if, in decreasing order of importance, the cluster 
virial mass at fixed X-ray temperature is being underesti- 
mated, S c is lower than the canonical value 1.7 and /, the 
assembled fraction of a cluster virial mass after which the 
X-ray temperature does not change significantly, is assumed 
greater than 0.75. However, the single most important fac- 
tor in determining the theoretically-expected overall value 
for N(> 6.2 keV, 0.32) is the present-day normalization for 
the dispersion of the density field, erg, which in turn re- 
sults from the observational value for the present density 
7V(> 6.2 keV, 0.05). 

Although we worked with all X-ray clusters that make 
up the dataset in Henry (1997), and even estimated the ef- 
fect of also considering the 5 clusters with lower X-ray fluxes 
present in the EMSS in the redshift bin from 0.3 to 0.4, in 
fact we only used the abundance of clusters with X-ray tem- 
peratures in excess of 6.2 keV to constrain fio. We mentioned 
some of the reasons for this choice in Section 3. Nevertheless, 
we decided to repeat the same calculations for a threshold 
X-ray temperature of 4.8 keV. This value also well represents 
the mean curve going through the observed cumulative X- 
ray temperature distribution function at both z = 0.05 and 
z = 0.32. 

The results regarding the best-fit value for fio, presented 



in Figure 4, are somewhat different from those we obtained 
when the threshold X-ray temperature was assumed to be 
6.2 keV. This is particularly true if the correction for the 
possibility of any of the 5 clusters with the lowest X-ray 
fluxes in the 0.3 < z < 0.4 EMSS sub-sample having X- 
ray temperatures in excess of 4.8 keV is included, as can 
be seen in Figure 5 for the open case. While the standard 
analysis without these 5 X-ray clusters prefers a value for 
fio between 0.4 to 0.5, when the correction for the scatter 
in the relation between the cluster X-ray temperature and 
luminosity is included, in the way described in subsection 
3.3, the preferred value for fio decreases to about 0.3. Now 
the fin = 1 hypothesis is excluded at more than the 95 
per cent confidence level, with or without the correction. At 
the 90 per cent confidence level, one finds that fin > 0.8 is 
excluded without the correction, being this limit lowered to 
0.7 when the correction is included. 

One can also estimate the joint probability of some fio 
value being excluded on the basis of the results relative to 
either one or both X-ray temperature thresholds. Assuming 
the data used in the calculations for the two thresholds is 
independent, the results then imply that the favoured value 
for fin is close to 0.55 (0.50 if the incompleteness correction 
is included) and the fio = 1 hypothesis is excluded at the 
99 per cent l evel. This agrees very well with the results of 
Henry ( 1997 ) and Eke et al. ( 1998 ), leading us to believe 
that the main difference between our analysis and theirs is 
our decision to draw our conclusions solely based on the ex- 
clusion levels obtained for the X-ray temperature threshold 
of 6.2 keV. 

A further potential problem one must consider when 
working with clusters whose observed X-ray temperature is 
as low as 4.8 keV is the possibility that the energy in the 
intracluster gas has increased as a result of (pre-)heating 
by supernovae and starbursts in the cluster galaxies. In fact 
this is the leading hypothesis (e.g. Navarro, Frenk & White 
1995; Markevitch 1998) put forward to explain the discrep- 
ancy between the observed slope of the X-ray temperature- 
luminosity relation, close to 0.3, and the expected val ue of 
0.5 if clusters evolve in a self-similar way (Kaiser 1986). 

Following Eke et al. ( 199S ) , we assume that in a cluster 
whose observed X-ray temperature is 4.8 keV, 17 per cent of 
its energy, that is 0.8 keV per intracluster gas particle, was 
due to (pre-)heating produced by processes occurring inside 
the cluster galaxies. This is approximately the amount of 
energy that gets injected into the int raclus ter gas particles 
in the simulation of Metzler & Evrard ( 1994 ) , where a galaxy 
cluster's X-ray temperature, which would otherwise be 5.6 
keV, increased to 6.4 keV. Note however that in the scheme 



proposed by Eke et al. ( 1998 ) a cluster this large would not 
be (pr e-)heated to the extent simulate d by M etzler & Evrard 
(1994), as in their proposal Eke et al. (1998) assume that the 
energy gained by each intracluster gas particle due to (pre- 
heating decreases as a galaxy cluster becomes larger, being 
close to zero for galaxy clusters with X-ray temperatures 
exceeding 6.2 keV. 

The above assumption means that the observed values 
for N(> 4.8 keV, z), when z = 0.05 and z = 0.32, should 
now be compared with the theoretically-expected values for 
7V(> 4.0 keV, z) at those redshifts. The resulting exclusion 
levels on the value of fio can be seen in Figure 5 for the 
open case. There is little difference compared to the results 



© 0000 RAS, MNRAS 000, 000-000 



Galaxy clusters at 0.3 < z < 0.4 and the value of Qq 11 



Open, > 4.8 keV 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



Figure 5. The absolute exclusion levels for different values of f2o 
only for the open case, when the threshold X-ray temperature 
of 4.8 keV is used. The full curve includes a correction (FC) for 
the possibility of the 5 clusters with lowest fluxes in the EMSS 
located between z = 0.3 and z = 0.4 having X-ray temperatures 
in excess of 4.8 keV. The dashed curve includes a correction (HC) 
for the possibility of (pre-)heating of the intracluster medium due 
to processes within the cluster galaxies. The dotted curve includes 
both corrections. 



in Figure 4 that follow from the standard no-heating cal- 
culation. The lower value for erg, required to match theory 
and observations at z = 0.05, more than compensates for 
the expected increase in the number of galaxy clusters with 
X-ray temperatures in excess of 4.8 keV at z — 0.05, in 
effect bringing this number down. In fact, the standard no- 
heating calculation for a threshold X-ray temperature of 4.8 
keV requires a value for erg (fio), so that the observed value 
for N(> 4.8 keV, 0.05) is reproduced, that is less than 3 per 
cent below that required by the > 6.2 keV data, quoted in 
equation (||) . On the other hand, including the (pre-)heating 
correction, the required erg (fio) value drops to 19 per cent 
below that preferred by the > 6.2 keV data. Though the 
coincidence between the erg values obtained for the two X- 
ray temperature thresholds 4.8 keV and 6.2 keV under the 
no-heating assumption may be accidental, it could indicate 
that (pre-)heating was relatively unimportant at least for 
the galaxy clusters observed at z = 0.05 with X-ray temper- 
atures exceeding 4 keV. If (pre-)heating was more impor- 
tant in the past than today, then the required erg (fio) value 
would be that obtained through the standard no-heating hy- 
pothesis, but the comparison at z = 0.32 would include the 
(pre-)heating correction. This would push the theoretically- 
expected value for N(> 4.8 keV, 0.32) up, favouring higher 
values for fio. This is not as far-fetched as it may seem, given 
that it is well known that the star-formation rate peaks be- 
fore z — 1 (e.g. Madau, Ferguson & Dickinson 1998; Baugh 
et al. 1998) , and consequently so does the rate of supernovae 
Type II (the rate of supernovae Type la peaks a few Gyr 
later) and the probability of starbursts. 

The results for the 4.8 keV threshol d X-r ay temperature 
are close to those found by Eke et al. ( 1998j ), leading us to 
believe that their exclusion levels for fio are dominated by 
the information associated with the threshold X-ray tem- 
peratures 4.0 keV and 5.0 keV. In our view, the analysis 
for these X-ray temperature thresholds carries with it a suf- 
ficient number of uncertainties, due to the problems men- 



tioned above, so as to render the constraints imposed on fio 
not very trustworthy. Only the data regarding clusters with 
X-ray temperatures in excess of about 6 keV seems suffi- 
ciently free of modelling problems so as to be potentially 
useful in constraining fio. 

Another possible complication has a risen from recent 
work by Blanchard, Bartlett and Sadat (1995) who use a 



sample of 50 galaxy clusters with mean redshift of 0.05, 
which were identified through the ROSAT satellite, to es- 
timate the cumulative X-ray temperature distribution func- 
tion at z = 0.05. They claim the number density of galaxy 
clusters at z — 0.05 with X-ray temperatures exceeding 4 
keV is being underestimated when the Henry & Arnaud clus- 
ter sample is used. Through the X-ray cumulative temper- 
ature distribution function at z = 0.05 they obtain, they 
then estimate fio using the EMSS cluster abundance in the 
redshift bin 0.3 < z < 0. 4 and the X-ray temperature data 
gathered in Henry (1997). They find the favoured value for 
fio to be 0.75, while fio < 0.3 is excluded at more than 
the 95 per cent level. These results coincide very well with 
ours when only the 6.2 keV threshold X-ray temperature 
is considered, thus perhaps implying that the discrepancy 
between the favoured value for fio found when different X- 
ray temperature thresholds are considered may arise from a 
underestimation of the cumulative distribution function at 
z = 0.05 for X-ray temperatures below about 6 keV. 

Unfortunately, due to uncertainties associated both 
with the observational measurements and the theoretical 
modelling of cluster evolution, the presently-available data 
on galaxy clusters with X-ray temperatures exceeding about 
6 keV is not able to strongly discriminate between cosmolo- 
gies with different values for fio. And in any case, the data 
available is probably not yet statistically significant. More 
is needed to support or disclaim the preliminary conclusions 
that can be obtained from it. In particular there are some 
oddities with the sub-sample of EMSS galaxy clusters ob- 
served by Henry, such as the strange redshift distribution, 
strongly clustered around 0.32, and the unexpectedly low 
X-ray temperature of MS2137.3, that makes one have some 
doubts about how representative this dataset is of the Uni- 
verse. 

Within the next few years, with the launch of the XMM 
satellite, possibly in late 1999, a significant increase in the 
quantit y and qualit y of the available data is expected to 
occur ( Romer 1998 ). It should then be possible to place 



stronger constraints on fio on the basis of the evolution of 
the galaxy cluster X-ray temperature function. This would 
be helped by improvements in the theoretical modelling of 
cluster evolution, perhaps based on the high-resolution hy- 
drodynamical iV-body simulations on cosmological scales ex- 
pected in the near future. 
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